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Abstract 



An exactly solvable model of two-component interacting Fermi vapour 
in two dimension within Thomas Fermi approach has been proposed. We 
assume a realistic off-diagonal s-wave interaction between fermions in the two 
hyperfine states. The interaction is taken to be given by a screened Coulomb 
interaction which is both relevant and is analytically tractable. There are two 
distinct limits in the case of trapped fermionic systems: When the rotating 
frequency is less than the trap frequency a Thomas-Fermi approach may be 
used for a reasonably large number of particles. In the case of rapidly rotating 
fermions, an alternative variational approach is used, where the local density 
approximation may be applied after projecting the system on to the lowest 
degenerate level. Analytic expressions for the relevant physical quantities are 
obtained and their significance is discussed in both these cases. 
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I. INTRODUCTION 



The experimental observation of Bose-Einstein Condensation(BEC) in extremely cold 
dilute gas of atoms has generated a lot of activity both on experimental as well as theoretical 
fronts on the properties of trapped dilute gases near the quantum degeneracy limit. There 
is now a renewed focus on the properties of trapped dilute gas of fermionic atoms at low 
temperatures. Magneto-optical confinement of fermionic gases has been reported for 6 Li |I[ 
and 40 K 0. A recent significant success in this direction is the experimental observation of 
quantum degeneracy in a dilute gas of trapped fermionic atoms [Q. They have realised the 
magneto-optical trapping of m K atoms in two different hyperfine states corresponding to 
\F = 9/2, mp = 9/2 > and \F = 9/2, nip = 7/2 > and finally reaching a single component 
Fermi vapour by selective removal of atoms in one of the hyperfine states. Amoruso et.al. 
H have studied the ground state and small amplitude excitations of such two component 
degenerate Fermi vapours placing them in spherically symmetric trap and discuss the effect 
of interaction on dissipative hydrodynamic processes in them. Salassnich et.al. jj] have 
studied the thermodynamic properties of such multicomponent fermi vapours theoretically 
and have pointed out that such systems can show phase separation and as well as shell 
effects. 

While the effective interaction between 40 K atoms in two different hyperfine states is 
repulsive a negative scattering length for 6 Li atoms holds promise of achieving a superfluid 
state in a mixture of 6 Li atoms prepared in two hyperfine states ||. This has motivated a 
number of theoretical papers. Recently collective excitations of the system in the normal 
phase [U and in the superfluid phase || have been investigated. Ghosh and Sinha [|l(J have 



studied |10[ the low lying collective excitations of Fermi and Bose systems in an anisotropic 
trap. 

Especially relevant to this paper is the recent interest in the properties of rotating bose 
condensates |Tl| and the experimental observation of vortices in stirred BEC of 87 Kb |12 



The properties of ground and low excited states of a rotating and weakly interacting Bose- 
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Einstein condensate in a harmonic trap has been investigated recently in various limits 
p^T5[. Benakli et.al. Jjlj] have studied the properties of macroscopic angular momentum 



state of a rotating BEC in a toroidal trap and identified /-dependent density profile as 
a signature of condensate rotation and superfluidity. Another very interesting theoretical 
issue related to the rotating BEC is the exact mapping of the problem in the limit of weak 
interaction to that of Lowest Landau Level problem in Quantum Hall effect |17| • Recently, 



experiments [0] has been able to trace a collective state of 160 vortices in a BEC consisting 
of Na atoms by imparting a large angular momentum to the system by rotating it. . This 
experiment is a major step towards realising this weakly interacting regime where the LLL 
mapping is useful. Theoretical model of such vortex lattice state has been discussed by 
Ho by invoking the fact that this system has a strong resemblance with the Quantum Hall 
systems. ( fl9|l ). 

It is therefore quite natural to extend the analysis to Fermi vapours also. Ho and Ciobanu 
fl8| have recently examined theoretically the nature the ground state of a rapidly rotating 
trapped Fermi gas in two and three dimensions in the noninteracting limit. They show that 
the density profile acquires features reflecting the underlying Landau level like energy spec- 
trum when the frequency of rotation approaches the bare confinement frequency. Properties 
of the spectrum of such a system had been earlier investigated in detail by Bhaduri et al 

A simple and straight forward way of studying the ground state properties of such in- 
teracting many-body system is through the Thomas-Fermi(TF) approximation. Butts and 
Rokshar studied the momentum and spatial distribution of the noninteracting system 



and identify the region in which TF approximation is valid. The ground state properties and 
the addition energy spectra of a two-dimensional interacting Fermi system has been studied 



by Sinha et al |22] within Thomas Fermi approximation using an exactly solvable model 
interaction. A similar exactly solvable TF model has been discussed for a two dimensional 
parabolic quantum dot in a weak magnetic field by R. Pino ( 

In this paper we analyse a multicomponent Fermi vapour confined in a rotating trap in 



two dimension within the Thomas Fermi approximation and with an effective interaction 
between the fermions. If the confinement along the axis of rotation (z- axis) is relatively 
much stronger, then only a few lowest energy subbands along that direction get populated. 
In such a situation the system becomes effectively a multicomponent two dimensional Fermi 
gas, the number of components being the number of subbands being occupied. In particular 
when one uses optical traps such a system is experimentally realisable [pip - In the case 



of rotating Bose gas such two-dimensional systems have been extensively studied Hl5| , but 
similar studies for their fermionic counterparts is still far from complete. The effective 
interaction is given by a screened Coulomb potential with a finite range. We show that 
within the TF approximation the system may be solved exactly and discuss the consequent 



results. We follow closely the formalism developed by Gallego and Das Gupta |25j who 
developed the TF method for the rotating nuclei. An attractive feature of the TF method 
is the ease with which non-trivial many-body solutions can be obtained even when they are 



interacting. As we have shown in a preceeding work [2^|, in some cases the results may be 
obtained purely by analytical methods. In this paper we discuss a model interaction which is 
relevant to physical systems. Appropriate limits can reproduce the results already obtained 
in ref. PEfl . Though our discussion is centered around rotating fermions in two dimension 
it is equally applicable to systems like quantum dot in a parabolic confinement with 



or 



without a magnetic field [22 



The paper is organised as follows. In section II we outline the derivation of the Thomas 
Fermi energy functional for a system of interacting rotating fermions in parabolic confine- 
ment and discuss our choice of effective interaction. In section III we apply the formalism 
to a two-component Fermi gas when the rotational frequency is less than the confinement 
frequency and obtain analytic expressions for various quantities of physical interest. We also 
present numerical results for the interaction energy, the rotational energy and the spatial 
density in a given angular momentum sector for a given number of particles. Furthermore 
we discuss the extreme quantum limit where a rapidly rotating Fermi gas can be mapped 
exactly to the problem of a charged particle in uniform transverse magnetic field in 2 - 



dimension. Exploiting this similarity we have analysed the many particle spectrum through 
a variational technique which is an improvement over the usual semiclassical method. The 
last section contains a discussion and summary of the results. 



II. THOMAS FERMI ENERGY FUNCTIONAL 

We derive the Thomas- Fermi (TF) energy functional for a confined two dimensional ro- 
tating fermi gas starting from the microscopic Hamiltonian. We shall follow the approach 
outlined by Gallego and Das gupta |25| who derived the TF density functional for rotating 



nuclei, but applied to a two dimensional system. 

The microscopic Hamiltonian of a two-dimensional rotating Fermi gas of N particles 
confined in a harmonic potential is given by, 

N D 2 1 



muj n r 



o' i 



i=l i<j 

-^mwfriXp;], (1) 

where V(fi,fj) denotes the two body interaction between fermions, uo and uj denote the 
confinement and rotational frequencies respectively. 

We want to minimise the energy of this system subject to the constraints 



dfp(r) = N, (2) 

J*wmp*t> = L, ( 3 ) 

where N and L are total particle number and the total Angular momentum, f(f,p) is the 
semiclassical phase space density which we shall define later. The configuration space density 
p is obtained by integrating the phase space density over the momentum space. 

Using Lagrange multiplier p and u for each constraint respectively and regrouping the 
momentum dependent term we obtain the following expression for the energy functional 

/p2 
drdpf(r,p)[— - u(xp y - yp x )} 

+ J dr\V{f) - p]p(f) + pN + uoL (4) 



where the mean-field one body interaction V(r) is 

V(r) = ^mu 2 r 2 + V H (r). (5) 

Note that the formalism itself is valid for arbitrary confinement potentials as long as the 
meanfield has a smooth behaviour. However in the later sections we restrict our analysis to 
parabolic confinement. 

The Hartree term is given by, 

V H (r) = J d 2 rV(r,r')p(r r ). (6) 

We ignore the exchange correction as it is in general subdominant compared to the Hartree 
term. 

The semiclassical phase space density is given by 

2 

where we have accounted for a factor of 2 due to spin degeneracy and 

p 2 

e(p, r) = ^ + V(r) - u(xp y - yp x ) (8) 

denotes the energy density in phase space. We may rewrite this by completing the square 

as, 

e(p, r) = [(p x + muy) 2 + (p y — mux) 2 ] + V(r) muj 2 r 2 . (9) 

Zi ill Zi 

The shifted momentum simply indicates that the center of the Fermi sphere at any given 
point r in the rotating frame is displaced from the usual p — 0. This shift is of no relevance 
within the classical TF approximation except for the appearance of the extra centrifugal 
term in the energy functional. 

The TF energy functional is therefore given by, 

r nh 2 o 2 1 

E[p] = / d 2 r [— i- + -ratlr 2 + V H - fi]p(r) +fiN + cuL (10) 
J 2m 2 



where 

n 2 = ul-u 2 (11) 

denotes the effective frequency in the presence of rotations. The TF equation for the spatial 
density is obtained by a variation of the energy functional, namely, 

l^R + i m n 2 r 2 + V H = fi (12) 
m 2 

To obtain the ground state density in TF approximation one has to solve the above equation 
self-consistently with the boundary condition that the density vanishes beyound the classical 
turning point ro, that is 

p(r) = 0, r > r . (13) 
Note that in the absence of rotations Q = uj$ and the solution of the above equation describes 



the ground state in the absence of rotations which has been analysed in detail before |22 
We may remark that in general the above equation is difficult to solve analytically except 
in specific cases which we shall discuss later. 

The angular momentum carried by the system is given by, 

L = jM P f(r, pl(m - yp.) = nuftntiy (14) 

which is simply the classical expression for the angular momentum. We may therefore 
identify the Lagranges multiplier uo with the angular frequency. 

We may also define the related energy functionals. The free energy of the system is given 

by, 

F = E — fiN — uoL (15) 

The TF equation fll2"|) may also be obtained just by varying the free energy. We also define 
the energy in the rotating frame through 

E = E' + uL, (16) 



where E is the total energy in the static frame. With the help of eq. fll2|) it can be shown 
that irrespective of the nature of the interaction E' is given by, 

E' = -fxN + ^mtt 2 J d 2 r r 2 p. (17) 

Note that it is necessary that we choose uo < ujq as otherwise the system becomes unstable. 

This concludes the basic TF formalism applied to rotating fermionic systems. In the 
following section we solve the TF equation for a specific case. In an earlier brief note ( p6| ) 
we considered some examples of model interactions for which the TF equations are exactly 
solvable for a circularly symmetric density profile. In here we extend our calculation to the 
case of effective (screened) Coulomb interaction. Depending on the range (which is a free 
parameter in our calculation) the interaction yields both logarithmic as well as the contact 
interaction. 

The Screened Coulomb interaction for a two dimensional system is given by 

where g is the coupling constant and the S denotes the relevant mass scale which determines 
the range of interaction. The choice of the effective interaction is dictated not only by the 
fact that it is tractable, but also by the fact that it is a good approximation to the realistic 
interaction between charged fermions in two dimensions. Given the very general form we 
have used in terms of the coupling constant g and the range b even for the neutral atomic 
(Fermi) vapours the consideration of such finite range interaction is quite important ( [[29|j). 
It can be checked that the above integral can easily written as 

V(r)=gK (~), (19) 

where K is the modified Bessel function and b = | denotes the range of the interaction. 

When the range is very small then the given form of V(r) behaves almost like contact 
interaction and in the limit of b — ► oo the interaction goes over to the logarithmic interaction 
which is Coulomb interaction in an ideal two-dimensional system. These can be checked from 
the following limiting behaviour of the modified bessel function K (z) |28| : 
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jr (£)~-bg(£) (20) 



and 

I -> oo 

- 27r6 2 5 2 (r-0 (21) 

Hence in the two extreme limits we obtain the logarithmic and the contact like behaviour. 
In Fig. 1 we have plotted for comaprisons three different potentials V(x) as a function oix 
- namely the Coulomb in two dimension 1/x, the negative of logarithmic interaction log\x\, 
and the screened Coulomb K (x). All distances are measured in units of length is defined 



by — y ^ • For the purpose of comparison we have set the range 6=1. It is easy to see 
that at long distance screened Coulomb falls faster than the unscreened Coulomb potential 
whereas at the sort distance it goes like -log\x\. We also note that for comparison we should 
only look at the portion of the plot of —log\x\ upto x = 1 where the potential vanishes. 
The TF analysis in these two limits has already been done in the previous paper |26||. In 



the following section we solve the TF equations for a general screened Coulomb interaction 
and discuss the results. The only approximation which has gone into obtaining the solutions 
is the assumption that the density is circularly symmetric. 



III. TWO-COMPONENT FERMI GAS WITH EFFECTIVE SCREENED 

COULOMB INTERACTION 

For a two component system the coupling constant can be written in a 2 x 2 matrix 
form to take care of inter and intra component interaction. We assume that such interaction 
strength matrix is symmetric in its off-diagonal components i.e., g\i = #21 — 9- Due 
to the anti-symmetry of the many body wave function the fermions occupying the same 
hyperfine state do not have the s-wave interaction but interact through the much weaker 
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p-wave interaction. To obtain the leading order results, we neglect this. Therefore we set 

011 = 922 = 0. 

The Thomas Fermi free energy functional of such two component Fermi vapour is there- 
fore given by 

F[ Pl ,p 2 } = I dr[—{p\ + pi) + imfiV( Pl + p 2 ) 
J m 1 

+ Pl (f) J dr'V 12 (\r- f\)p 2 (f)} - HiNi - p 2 N 2 , (22) 

where p±, p 2 denotes the densities of the two components while Ni and N 2 denote the pop- 
ulations of the two hyperfme states which is kept fixed during minimisation. Furthermore, 
we have 

_ g 12 f e ik -F-^dk 
V ^ r - rl)= 2^J -J^T5^- (23) 

for the interaction matrix. 

In what follows we use the dimensionless variables (denoted as variables with ~ on top). 
The quantities with energy dimensions are measured in units of huo and the length is 
measured in units of l = \J After this rescaling the free energy functional may be 
written as, 

F[p ± , p 2 \ = J dx[n(pl + pi) + ^ V(P 2 1 + P2) + ~ Pl J dx'V 12 (x, x')p 2 (x')} - p.N, - p 2 N 2 , 

(24) 

where scaled variables are given by 

p a = 



^0 

Pa — Pa^o 

X = y (25) 

Minimising the free energy functional, we get the equations for the TF densities, 
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where 



2u~ Pl = fi 1 - V m - ^Q 2 x 2 , (26) 



2vrp 2 =fa- V H2 - -OV, (27) 



V Ha {x) = j^J dx'V aP {\x - x'\)p p {x') (28) 



By definition the Hartree potential obeys the following set of Poisson's equations, 



V 2 V Hl = -2ngp 2 + ^M, (29) 



and 



V 2 V H2 = -2ir~g~Pi + (30) 

b 2 

where the coupling g = ^ and the range parameter is now given by b = ^ and, pi j2 is the 
density in dimensionless units (p = p(r)lo) 



Substituing Eqs. (|29|j30|) in eqs. (|26| , |27|) we get following equations for the two densities, 



2vrV 2 pi + 2tt 2 - 2ngp 2 + ^(fii - 2np 1 - \x 2 & 2 ) = (31) 

b 2 2 

and 

2nV 2 p 2 + 2tt 2 - 2ngp 1 + i(p 2 - 2np 2 - -x 2 Q 2 ) = 0. (32) 

b 2 2 

The eqs. ( pi| ) and fl52|) may be solved subject to the usual boundary conditions 

pi = at x > x\ (33) 
p 2 = at x > x 2 (34) 

At this stage it is useful to define the following quantities 

P = Pi+P2, (35) 

s = p x - p 2 , (36) 

= (pi ± p 2 ), (37) 
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where p denotes the total density, and s denotes the spin density. It is easy to see that in 
terms of the density and the spin density the coupled equations for the two components are 
diagonal. In this paper we consider the case when the system is unpolarised, that is when 
the spin density is zero. The case of partially polarised system will be discussed in detail in 
a forthcoming publication |27J . 

Before proceeding further we note that, since Q 2 = uj 2 — uj 2 , the system becomes unstable 
when the rotational frequency approaches the confinement frequency. TF method can be 
applied abinitio only to the case when the rotational frequency is less than the confinement 
frequency. In the rapidly rotating case many levels collapse into the lowest one reminding 
one of the Lowest Landau Level(LLL) structure. We discuss these two cases separately in 
what follows. 



A. Trapped rotating fermions 

When the particle number in the two hyperfine states are equal the system is unpolarized. 
Since the interaction strength is symmetric, the turning points are therefore the same. This 
situation changes only if there is a spontaneously broken symmetry in the ground state 
leading towards phase separation. 

For the unpolarised case the spin density s(x) vanishes everywhere. Adding eqs.flST]) and 
([321) we obtain 

2vrV 2 p + AQ 2 - 2irgp + i(/I + - 2vrp - x 2 Q 2 ) = (38) 

b 2 

The same equation may also be obtained by minimising the energy functional given in the 
previous section, see eq.(|T0P where we simply considered a system with spin degeneracy of 
2 as in the unpolarised case. 

The solution of the above equation may be written in terms of the modified Bessel 
function, namely 

2np = AI {rix) + ^ + + - x 2 n% (39) 
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where 

v 2 = (g + i) (40) 

b l 

The solution for the density has three unknown parameters, the turning point xq, the overall 
normalisation A and the chemical potential ft + . In the following we outline how these may 
be determined from the known constraints. The first constraint arises from the fact that at 
the classical turning point the density vanishes. We may use this condition to determine the 
chemical potential: 



A + = -bWAIo(vxo) + d 2 xl - (41) 



This gives /i + in terms of the other two unknowns A and x . Substituing this we get 

^- AI ^M)- l]+ w [A - A (42) 

Interestingly, the first term in the above expression corresponds to the density obtained in 
the case of logarithmic interaction and the second term corresponds to the density obtained 
in the case of contact interaction if the coefficients are properly adjusted. These different 
limits were derived earlier in Ref. separately for these two interactions. 

The second constraint is given by the continuity of the Hatree potential and its derivative 
which is the electric field experienced by the charge distribution at the turning point. Adding 
eq. (|29|) and eq. (0) we have 

V 2 V^ + = -2^p+^|±, (43) 
b 2 

where Vh+ corresponds to the Hatree potential with 

V H + = Vhi + V H 2 

yln the region x > xq the equation ( |4*3| ) takes the form 

V 2 VW = %t (44) 
b 2 
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since the density is zero beyond the turning point. The solution is given by 

y H+ = W), (45) 

where A\ is an overall normalisation which is not needed in what follows. The electric field 
is just the negative of the gradient of this potential E = — VV. We define the ratio of the 
potential to the field as follows 

= bKoixoM/K^xo/b') (46) 

From eq. ( |26|) and (^7j) we have for x < Xq 

V H+ = fl + - Q 2 x 2 - p (47) 

Again taking the ratio of the potential and the electric field we have 

V H +(x ) _ p - tt 2 xl - p(x ) 



E(x ) 2fi 2 x + {p{x))'\ x=X0 (48) 

Now demanding the continuity of the potential and the field across the turning point we 
have from eqs. fl46|) and fl4~8|) 



2n 2 fb 2 (2 + 



xoKo(xo/b) ■ 
bK^xo/b) 



ai ( V x ) = --^f- — : — t; m > ( 49 ) 

v> v +vnl Ioivxo)Kl( ^ ) \) 

which now gives the relation between A and the turning point. Finally the turning point, 
and hence all the undetermined parameters are determined using the fact that the number 
of particles in the system is fixed. That is, 

t]I (t]x ) 2 4 b 2 r] 2 

Note that the turning point increases with N in a nonlinear fashion given by the above 
relation. For large number of particles, however, the turning point is given by the following 
asymptotic relation 

1 2 r 4 

N^-Z^ (51) 
4 b 2 rj 2 v ) 
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All the other relevant quantities defining the system may now be written down. The exact 
expression for the total angular momentum is given by 

h lu r] 7] 2 I {rjx ) T] 2 4 uj b 2 r) 2 12 

The Energy is again given by the expression 

E + 1^1 



/ic^o 2 4 cucjo ft 
It is easy see that all the above expressions above can be written as a sum of contributions 
coming from the purely contact interaction and the purely logarithmic term ( [^6j), two 
extreme limits of the screened Coulomb interactions. However this is not valid for the 
interaction energy since it is bilinear in the density. This may be seen from the analytic 
expression for the interaction energy 



Ei n t E 2 f ~2 i 1 Q 2uJ L 

7T 



J p 2 xdx — -- 



huo huo J 2 u h 

= ^-4 (/l + /2 + /3) -2^' (54) 



where 



/, = 1{ AI M ^ - - 2*.£g + S], (55) 



J 2 = A/o r/x )~— 2— 1 — — 56 

b 2 ?] 2 V x r)I (r]x ) 8 

1 O r 6 

* ^ J'&T < 67 > 

While the system is analytically solvable as shown above, the actual behaviour is not 
often easy to see. We illustrate the results for reasonably large systems with number of 
particles iV = 50, 100, 1000. The numerical calculation is done in two different scenarios. 
In the first case we have kept the particle number N, the range b and coupling strength g 
fixed and vary the external frequency of rotation. In the second case we have calculated the 
relevant quantities as a function of the quantity gb 2 where all the other quantities except 
the range b are kept fixed. 
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In Fig. 2 we show the turning point function of the frequency to measured in units 

of ujq. Note that Eq.(|50D is highly nonlinear. In order to estimate the turning point we use the 
value given by the asymptotic result as an upperbound and use this value to iterate the exact 
equation to obtain the exact value of xq. It is easy to see that the system size increases with 
rotational frequency. All other quantities are determined once the turning point is known. 
The system size increases with rotational frequency rather slowly as expected. The size 
variation with respect to the number of particles is approximately given by the asymptotic 
formula. 

The angular momentum per particle also follows a similar behaviour, as shown in Fig. 
3. The ratio of the energy of the system measured in terms of the ground state energy 
is shown in Fig. 4. Shown also in this figure is the interaction energy for N = 1000. At 
low frequencies the interaction contributes substantially to the total energy. However, at 
large frequencies the interaction contributes little, and most of the total energy is simply 
the energy of rotation which is proportional to the angular momentum L. Following the 



analogy with rotating nuclei, this is usually referred to as the Yrast region 30J and has been 
a subject of intense study in recent times [T3[ . 



As already pointed out, the screened Coulomb interaction can mimick different types 
of two body interaction depending on the range. We study the system for N = 1000 at a 
fixed frequency of uj = 0.25a>o, in Figs. 5,6,7 where we have plotted the turning point, the 
angular momentum per particle and the total energy respectively as a function of gb. The 
quantity gb 2 here is analogous to the coupling strength in the pseudopotential of a three 
dimensional Bose condensate, g^n = 4n ^ a , where a is the s-wave scattering length. Note 
that the system of fermions, for these purposes, may be regarded as a hardcore bosons. This 
three dimensional coupling constant is related to the two dimensional coupling constant by 



the formula ( |T5| ) g 2 D = dsDy 111 ^- We have estimated this constant using the parameters 
for 40 K and chosen the appropriate range of gb 2 (see ref. [ID] for such an estimate). From 



these figures, it is easy to see that within the Thomas Fermi approximation the effect of 
change of range is considerable even at moderate rotational frequency as above. In particular 
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in Fig. 7, we see that both the total and interaction energy increase with the range. 

It should be mentioned that the analysis can also be extended to the partially polarised 
system where the number of particles in different components are not same . Formally one 
can write the solutions in closed form which however are more difficult to evaluate since the 
turning point for each component density will be different. This case will be considered in 



detail in another paper 27 



B. Rapidly Rotating fermions 

The above TF analysis is valid when the trapping frequency is larger than the rotational 
frequency. The effect of rotations is to reduce the effective trapping frequency and in the 
absence of interactions this leads to bunching of single particle levels. There are thus many 
similarities between the system of rotating fermions and the fermions in the presence of an 
external magnetic field along z-direction. The rotation frequency to plays a role analogous 
to the cyclotron frequency. We exploit this similarity in what follows. 

Noninteracting Hamiltonian for a rapidly rotating system is given by, 

ff o = E^ + ^ 2 -^. (58) 

i 

where as before uj is the external rotation frequency. The single particle energy levels are 
given by, 

E nJ = (2n + |/| + l)fnu - fuvl, (59) 

where n is the radial quantum number. This spectra is the same as that of a cranked 
harmonic oscillator with two frequencies ujq + uj and u)$ — uj (for details see Ref. [0). The 
Landau level like structure starts to form when (uj — uj) << (uj + uj) and the mixing 
between the different levels is supressed. In particular when uj — uj = 0, the spectrum 
becomes identical to that of a set of highly degenerate discrete Landau levels. In the other 
limit when the external rotation frequency is not too high and [i» (ujq + uj), so that many 
levels are occupied, the system is described by Thomas-Fermi approximation. 
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_ l2 9(R-r), (63) 



In the lowest Landau level (LLL) the particles have the radial quantum number n = 
and the angular momentum I > 0. The single particle wave functions and energies in LLL 
are given by, 

ifo = -^=(r Ye- r2/2l °, (60) 
Ei = (u — uj)hl + Tiujqi (61) 

where If. = and z = x + iy. 

The density of the non-interacting spin polarised electrons is given by, 

1 N 1 

Po = -7iEn(^o) i e- r2/ 'o, (62) 

7T£ l=Q l\ 
1 

where, R is the radius of the system for large number of particles. In this limit the particles 
form a flat droplet like state. 

We now consider the interacting two component unpolarised fermions in the "lowest 
Landau band" when the rotation frequency is almost equal to the external trap frequency. 
The interaction is taken to be the screened Coulomb interaction defined earlier. Since TF 
method can not be applied directly, we use the standard variational method for determining 
the ground state properties. We construct a local density functional by taking the expecta- 
tion value of the quantum mechanical many particle Hamiltonian with the variational state 
given above, neglecting the mixing between the ground and excited states. The beauty of 
this method is that all quantities of physical interest may be derived analytically. Though 
this method also does not goes beyond the mean field approximation , it performs better 
compared to the bare TF method since the variational ansatz does account for quantum 
mechanical effects. 

The many body variational wave function for each component may be written as, 

MM) = fazifte-Xi**™ (64) 

Kinetic energy of the system is then given by, 
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< T >=< V| - ^( 4 J^)I^ >= \rnulj r 2 p{r)d 2 r. (65) 

Note that the kinetic energy density is linear in p, that is, T[p] = ^muj 2 r 2 p(r) , where as in 
the standard TF case it is quadratic. The total angular momentum of the two-dimensional 
system is given by, 

/dip dip* f f 

[ip*z-— — zip] — ft I ip*ip = / d 2 r[mujQr 2 — h]p(r). (66) 
CJ As \J As tJ J 

Using these, the energy density functional of two component fermi gas can be written 

as, 

E[ Pl ,p 2 ] = T[ Pl ) + T[p 2 ) + V[ Pl ) + V[p 2 ] - Q(L[ Pl ] + L[p 2 }) + E mt [ Pl ,p 2 ]. (67) 
This energy functional is similar to the magnetic field Thomas Fermi method developed by 



Lieb et al |33[]. For the unpolarised system, p\ — p 2 . The free energy of the system may be 
written interms of p = p\ + p 2 alone. Using the expression for kinetic energy density and 
scaling to dimensionless variables as before, the free energy of the two component system is 
given by 

' - (1 — Co) [ d 2 xx 2 p+ \ [ d 2 xid 2 x 2 p(xi)Vi 2 (\xi - x 2 \)p(x 2 ) — (fi — u) [ d 2 xp, (68) 



hto J '4 



where to = — . 



Minimising the free energy for a fixed number of particles gives the following integral 
equation for the density: 

(1 - Co)x 2 + - J d 2 x'V l2 {\x - x'\)p{x') = (ft-Co), (69) 

As noted before, the screened Coulomb potential satisfies, the following differential equation, 

V 2 V H = -2K~g~p+y§. (70) 
tr 

While the equation system looks complicated, the solution for the circularly symmetric 
case is rather simple and is given by, 

19 



p=-(/i c -x 2 ), (71) 

71 

where 

1 — Co . . 

/i c = ^ + 46 2 . (73) 

One can easily recognise that the solution bears a resemblence to that of the lowest Landau 
level droplet for the rapidly rotating fermions. 

The size of the droplet (R) may be obtained by normalising the above density to yield 
the total number of particles N which is fixed. We have 

li c = R 2 

R = (2N/a) 1/ \ (74) 

Now we are in a position to compute all quantities of interest relevant to the rapidly 
rotating two-component system. Angular momentum of the system is given by 

L = = -a(2N/af 2 - N. (75) 

lb 6 

For a given angular momentum this equation determines the parameter a interms of the 
fixed parameters L and N. That is, 



2 N 3 
a — —- 



9(L + N) 2 ' 
The interaction energy is then given by, 



(76) 



2 M 3 nh 2 

E mt (L, N, g, b 2 ) = - - g AT \L + N(l - 2b 2 )]. (77) 



9 (L + N) 2 

The total energy in the rotating frame is given by, 



E(L, N, g, b 2 ) = (l-u) J d 2 xx 2 p + Co J d 2 xp + E int 



N + ljTTW [2L + N{1 - 2i2)] (78) 
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Variation of the total energy with external rotation frequency is given by 

E(u, N, g, b 2 ) = N + ^]~gb 2 (l-Cj)(2Nfl 2 - (1 - u)N(l + 2b 2 ). (79) 

Since a oc (1 — u), the angular momentum increases as the rotational frequency approaches 
the confinement frequency. In fact when oj = uj , the total angular momentum blows up, 
and the system resembles the LLL droplet. The total energy is given entirely by the energy 
of non-interacting LLL system since the interaction goes to zero. Total energy goes to the 
zero point energy as ~ y/1 — oj. 

IV. SUMMARY 

We have discussed a Thomas-Fermi model for a rotating two-component interacting 
Fermi vapour in two dimensions. We believe that the interaction, the screened Coulomb, 
is realistic with interesting limits as the range is varied. The results correponding to short 
range and logarithmic interactions may be obtained with suitable approximations. More 
importantly the system is analytically solvable within the TF mean field approach. We 
have neglected the diagonal component of the interaction including the dominant s-wave off 
diagonal interaction. We have first considered the case when the trap frequency is larger 
than the rotational frequency since one can use the TF method ab initio. However when the 
rotational frequency is much larger than the trap frequency we first projected the system on 
to the lowest degenerate level and used a TF like approach to derive the analytical result. 
It should be noted that, in general, it is rather hard and near impossible to obtain accurate 
results for large number of particles using exact diagonalisation. In contrast, the TF method 
is solvable in a realistic case as shown and provides a reasonably good approximation. An 
important drawback in the present analysis is that we have not included the exchange 
corrections. This however is obviated if the number of particles in the system is large and 
sufficiently dilute. 

To put the utility of the method in perspective, unlike in the case of superfluid bosons, 
the system of fermions has a rigid body rotation. A direct measurement of the angular 
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momentum in a system of bosons is very difficult. But in the case of fermions total angular 
momentum can be estimated just from a measurement of the moment of inertia. It is 
valid for low rotation frequencies, where the usual TF method is good. The radius of the 
cloud can be measured easily by observing the diffusion of the gas (by switching off the 
trap potential). In the rapidly rotating case, however, the total angular momentum is not 
proportional to < r 2 > due to the quantum effects being large and a direct measurement of 
the angular momentum would be difficult. The system is however more interesting because 
of its similarity to the LLL problem. It would be interesting to see if the system exhibits 
quasiparticle (quasi hole )like excitations which is similar to the vortex excitations in the 
corresponding Bose systems. 
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FIGURES 



FIG. 1. The spatial variation of three type of potentials namely ^ (dashed), — log\x\ (dotted) 
and the screened Coulomb Kq(x) (solid line) (where we have set b = 1) is shown. 

FIG. 2. The plot of turning point of the rotating fermions xq as a function of the external 



frequency lo. The unit of lo is the confinement frequency loq and the unit of length is 1$ = y 
The three curves correspond to N = 50 (solid line),^ = 100 (short-dashed) and N = 1000 (long 
-dashed). For all the curves plotted we have chosen b = 0.15 and g = 0.1. 

FIG. 3. The plot of angular momentum per particle of the rotating fermions L/N as a function 
of the external frequency lo. The unit of lo is the confinement frequency loq and the unit of L is 
ft. The three curves correspond to N = 50 (solid line),iV = 100 (short-dashed) and N = 1000 (long 
-dashed). Again for all the curves plotted b = 0.15 and g = 0.1. 

FIG. 4. The plot of energy of the rotating fermions as a function of external frequency cj.The 
unit of lo is same as that in Fig.3. The energies are normalised to the lo = case. The upper curve 
with dashed line shows the variation of the total energy and the lower curve (solid line) shows the 
variation of the interaction energy. The number of particles is N = 1000 and b = 0.15, g = 0.1 

FIG. 5. The plot of the turning point of the rotating fermions as a function of gb when b is 
varied. The number of particles N, g and the frequency to (in units of luo) are respectively given 
by 1000,0.1,0.25. 

FIG. 6. The plot of angular momentum per particle L/N as a function of gb 2 when b is varied. 
The unit of angular momentum and length is again fi and Iq respectively. The number of particles 
N, g and the frequency u (in units of 000) are respectively set at 1000,0.1,0.25. 

FIG. 7. The plot of energy as a function of gb 2 when b is varied. The continuous curve cor- 
responds to the total energy where the dashed curve corresponds to the interaction energy alone. 
The number of particles N, g and the frequency lo (in units of ojq) are given by 1000, 0.1, 0.25 
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